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Abstract. Orbital information for Aryabhata was computed and made available for 
telemetry and tracking schedules and for mission analysis. Taking an available state- 
vector once in a week from USSR, the orbital computations were made for the first 
four months using a numerical integration procedure. During the later period of the 
mission, an analytical method was adopted to reduce the computer time. Aryabhata, 
in a near circular orbit at 600 km and having an inclination of 50.68°, experienced 
a constant atmospheric drag due to which the semi-maior axis decayed by about 
0-5 km in one year. The nodal regression rate, due to asphericity of the earth, was 
almost constant at 4-632°/day. The inclination remained constant throughout the year. 
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1. Introduction 

Aryabhata, the first Indian satellite, was launched into a 600 km near-circular orbit 
at 50-68° inclination on 19 April 1975 by a Soviet rocket from Cosmodrome. Orbital 
information for Aryabhata was computed using the programs developed at ISRO 
Satellite Centre. The main objective of the orbit computations was to provide 
data for telemetry and tracking schedules and for mission analysis. 

For the first four months, the satellite position and velocity vectors were computed 
accurately, by using a numerical procedure, based on an initial state-vector available 
(once in a week) from USSR. An analytical approach, for which prediction accuracy 
was lower but sufficient for telemetry and tracking schedules, was adopted later to 
reduce the computer time. Also the state-vectors were not available on a continuous 
basis after the first four months. 

The limited tracking data available from SHAR could not be used for the refine¬ 
ment of orbital parameters as the orbit determination programme was not ready 
for operational use at that time. However, the Doppler data were used to adjust 
the absolute time (equatorial crossing time) and occasionally the right ascension of 
ascending node £2, during the later period when analytical procedure was adopted. 

In both the programs only the effects due to the asphericity of the earth and at¬ 
mospheric drag are considered. It is assumed that all other perturbing forces cause 
negligible effects on the satellite’s motion. A brief description of the programmes 
and their performances is given here. 


2. Orbital computations 

The various external forces acting on a satellite may be classified into two groups, 
♦A list of symbols appears at the end of the paper 
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namely the central force and the perturbative forces. The central force is the gravi¬ 
tational force of a spherical earth which causes the spacecraft to move around the 
earth in an elliptical path. The other forces like the atmospheric drag, asphericity 
of the earth, luni-solar gravitational attraction, solar radiation pressure etc., consti¬ 
tute the perturbative forces. They are quite small in magnitude as compared to 
central force. These forces influence the path of the satellite (in an elliptic orbit) 
which is no longer a closed one (Escobal 1965; Herrick 1971). For close-earth 
satellites, luni-solar perturbations and solar radiation pressure (leaving some special 
cases) effects are several magnitudes less than those due to the asphericity of the 
earth (Robertson 1968). In the analysis of the orbit of Aryabhata, only the effects 
due to the asphericity of the earth and the atmospheric drag are considered. 

2.1. Numerical method 

The differential equations of motion with respect to a geocentric rotating co-ordinate 
system (x pointing towards the Greenwich meridian, y perpendicular to x axis in the 
equatorial plane and z axis pointing towards the north pole) including perturbing 
forces are (EPyasberg e t a l 1974) 


X = — - X + 2(0^ + ( 0 *X — Spvx + AiX, 

r 

■y = — P y — 2(o t x -f w*y — spvy + Ai>\ 
r 

z = — - z — spvz + A x'z, (1) 

where s=C D A/2m; spvx, spvy and spvz are the perturbing acceleration components 
due to the atmospheric drag; and Ai*, Aiy and Ai'z are the perturbing acceleration 
components due to the asphericity of the earth. 

The geopotential of the earth is given by 

U = p/r + U d , (2) 

where U d is the disturbing potential given by 

^ = r [ “ 2 = 2 ( tT 7 " P " (SiD + ( t )^ 2 - 2 (siD ^ 

cos 2A -f S 2t2 sin 2A^ J. (3) 

An analytical model proposed by EPyasberg et al (1974) was used to compute the 
atmospheric density, given by 

P=P H K 1 K 2 K 3 K i , 


( 4 ) 
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The differential equations of motion were integrated numerically by a 4th order 
Adams-Moulton predictor-corrector method to which the initial four values are 
supplied through a Runge-Kutta-Gill routine. Normally, computations of position 
and velocity vectors for a week with optimum step-size (40 s), took about 
1-5 hr of computer time on an IBM 360/44. The optimum step-size was chosen 
after studying the behaviour of the time period with different step-sizes for about 50 
revolutions. The behaviour is shown in figure 1. At the end of the week the predict¬ 
ed state-vector was compared with the current state-vector supplied by USSR. The 
difference was found to be of the order of 1 to 3 km in position and 1-3 s in the 
absolute time. One such example is presented in table 1. 

2.2. Analytical method 

As a first approximation for the satellite’s motion the Keplerian orbit is assumed. 
The orbit is represented by six parameters viz. the semi-major axis a, the eccentricity 



Figure 1 . Behaviour of the time period for 50 revolutions with different step-sizes. 




Table 1 . Orbit prediction programme 
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e, the inclination i, the argument of perigee <o, the right ascension of ascending node 
£2 and the equatorial crossing time T. In due course, the variations in the orbital 
parameters due to the perturbing forces are added. Only the secular effects due to 
second zonal harmonic (/ 2 ) and the atmospheric drag are considered. No account 
is taken of long and short periodic perturbations. The expressions for the varia¬ 
tions of orbital parameters are taken as given by King-Hele (1968), Sterne (1968) 
and Robertson (1968). Due to the many approximations made in the orbit model, 
the accuracy of the predictions was lower. At the end of the week, when the 
predicted state-vector was compared with the current state-vector supplied by 
USSR, the difference was found to be of the order of 50-100 km in position and 
2-3 min in the absolute time. However, the accuracy was sufficient to keep the 
telemetry and tracking schedules. The important advantage of this method is that 
it takes only 10 min of IBM 360/44 time for predictions over a week’s period. 

2.3. Computed orbital information 

2.3a. Target designation 

The information from onboard the satellite can be extracted by the ground station 
when the satellite is above its horizon. The time spent by the satellite from horizon 
to horizon is called visibility period. When h is greater than 0°, the satellite is 
above the horizon of the station considered (figure 2). A control function to 
check the visibility condition is 

CF-(^-sta*4 (5, 


Figure 2. Geometrical representation of a 
satellite in a circular orbit around the earth 
for calculating the visibility period. 
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The satellite rises above the horizon of the station when CF changes sign from 
negative to positive, and sets when sign change is from positive to negative. Thus 
the visibility is determined. 

During the visibility period, the position of the satellite in the topocentric 
co-ordinate system was determined by a co-ordinate transformation. The origin 
of the system is the ground station, and the axes are towards the local north, east 
and zenith directions. The position of the satellite with respect to the station is 
given by 


P T = i — G. (6) 

The co-ordinate transformation is accomplished as follows. 


Pn 


sin <j> j cos 6 

—sin <f> s sin 6 

cos <f> s 


~ Px 


Pt 

= 

—sin 6 

cos 6 

0 


Py 

(7) 

_Pz_ 


cos <f>, cos 6 

cos <j> s sin d 

sin& _ 


-P*- 



where 0 = 0 g + A s . 

From the above transformation the azimuth and elevation are evaluated by 
A, = tan" 1 ( P Jp n ), h = sur 1 (pJp T . 

The range, azimuth and elevation history at 30 s interval during the visible portion 
of the orbit are called target designation and are used for orienting the tracking 
antenna in the direction of the satellite. 

Aryabhata is visible to SHAR ground station four times a day. The maximum 
visibility period is about 13 min and the average visibility duration over a day is 
about 40-41 min. 

2.3b. Eclipse characteristics 

Since the power generated by the solar cells and the thermal condition of the onboard 
systems are dependent on the length of time the spacecraft spends in the sunlight, it 
is important to know this time. The time spent by the spacecraft in the earth’s 
shadow can be computed by knowing the position of the satellite and the sun. Stan¬ 
dard methods are available for determining the sun’s position. The earth shadow 
is assumed to be cylindrical as shown in figure 3. Depending on the position of the 
sun, the shadow period for Aryabhata varies from 0 to 37 % of an orbit. For the 
first six months Aryabhata spent about 300 hr in sunlight. The eclipse characteristics 
of Aryabhata are shown in figure 4. 

2.3c. Ground trace 

Subsatellite point of a satellite is the pierce point of a ray dropped from the satellite 
in such a direction that it intersects normally the surface of the earth. The locus of 
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Figure 3. Geometrical representation of a satellite in a circular orbit showing 
the region of shadow. 



all the sub-satellite points is called the ground trace. Knowledge of ground trace is 
essential to know the region of the earth on which the satellite passes at a particular 
time. This information along with the height (height of the satellite above mean 
geoid) is used in attitude determination and to correlate with the observed data of 
the experiments. The history of the sub-satellite point is given by geodetic latitude 
(</>) and geographic longitude (A). The ground traces of Aryabhata for the first few 
revolutions are shown in figure 5. 


3. Tracking 

The tracking data of Aryabhata from Doppler, tone range and interferometry systems 
were available from SHAR. These data were smoothed by a pre-processing 
program (a preliminary version) developed at ISAC. These data could not be 
utilised for the refinement of orbital parameters as the orbit determination program 
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Figure 5. Ground trace of Aryabhata for the first 15 orbits after the launch. 

was not ready for operational use at that time. However, the Doppler data were 
used to adjust the absolute time and occasionally the right ascension of ascending 
node O. The observed minimum range, the calculated minimum range and the 
corresponding timings were compared and the differences were noted. The differ¬ 
ence in time is used to adjust the equatorial crossing time of the corresponding orbit. 
If too much difference is observed in the numerical values of the minimum range, 
the right ascension of node £2 is adjusted by a geometrical procedure. 


4. Performance of orbital predictions 

The main objective of the orbital computations for Aryabhata was to cater to the 
needs of telemetry and tracking schedules for which the requirements of accuracy 
were relatively less stringent. In the initial phase of operations (first 4 months) 
the orbital predictions were generated using the state-vectors available from USSR. 
The predictions were accurate. However, for generating the complete orbital infor¬ 
mation the program used to take about 9-10 hr of computer time on an IBM 
360/44 for a week. After 3 months of operation when the ground-trace information 
was suppressed, the program used to take about 3-4 hr of computer time. For 
this reason the analytical model (with many approximations) was adopted later, 
even though the accuracy of orbit prediction was lower. Also the state-vector was 
not available on a continuous basis from USSR for the later period. Even the few 
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Table 2 Orbital elements of Aryabhata 


Orbit No. 

Semi-major 

axis 

Eccentricity 

e 

Inclination 
i (deg) 

Right ascension 
of Node (deg) 

Argument 
of perigee 

W (deg) 

101 

6967-781 

0-003150 

50-707 

233018 

39-083 

201 

6967-744 

0 002918 

50-686 

191-849 

84-702 

301 

6967-611 

0-002730 

50-684 

160-791 

113-587 

401 

6967-730 

0-003087 

50-693 

129-828 

131-769 

501 

6967-735 

0-003057 

50-682 

98-787 

148-358 

601 

6967-761 

0-003099 

50-678 

67-790 

171-554 

701 

6967-673 

0-003235 

50-674 

36-779 

197-995 

801 

6967-757 

0-002820 

50-681 

5-774 

221-239 

901 

6967-599 

0-003481 

50-671 

334-765 

241-277 

1001 

6967-636 

0-003482 

50-689 

303-737 

265-973 

1101 

6967-632 

0-002921 

50-683 

272-807 

289-444 

1201 

6967-664 

0-002991 

50-686 

241-720 

318-314 

1301 

6967-742 

0-002820 

50-680 

210-703 

345-444 

1401 

6967-645 

0-003131 

50-678 

179-687 

16196 

1501 

6967-705 

0-002730 

50-683 

148-689 

40-430 

1601 

6967-725 

0-002881 

50-676 

117-680 

66-205 

1701 

6967-662 

0-002876 

50-685 

86-654 

90-996 

1801 

6967-695 

0 002932 

50-688 

55-309 

120-286 

1901 

5967-624 

0-02943 

50-686 

24-291 

145-030 

2001 

6967-531 

0 003948 

50-688 

355-322 

173-690 

2101 

6967-556 

0-003943 

50-688 

323-278 

198-249 

2201 

6967-493 

0-003955 

50-686 

292-253 

222-779 

3500 

6867-473 

0-003395 

50-695 

248-454 

173160 

5080 

6967-353 

0003427 

50-688 

118-782 

196-135 

5150 

6967-327 

0-003657 

50-688 

97060 

221-748 

5500 

6967-321 

0002286 

50.689 

348-540 

305-452 

5600 

6967-254 

0-002849 

50-687 

317.508 

329.887 


state-vectors available were received late. However, the main objective of the 
mission, viz. to keep the telemetry and tracking schedules, was successfully met. 
It was not possible to carry out a detailed orbital analysis for Aryabhata but 
the available state-vectors have been transformed to mean elements (short periodic 
perturbations subtracted) over a period of about one year, starting from orbit No. 
101 to orbit No. 5600, to know the history of the orbit. Very little information is 
available after the 2000th orbit as seen from table 2. Hence, this result should be 
considered as tentative. The following observations can be made regarding the 
behaviour of the different orbital elements. 

(i) The variation in the semi-major axis a over a period of about one year is shown 
in figure 6. The data points are approximated by a linear relation with time, show¬ 
ing that the semi-major axis for Aryabhata has decreased by about 0.5 km over one 
year duration. When the effect was evaluated by using a Jachia-65 atmospheric 
density model, the decay in a was found to be of the order of 1-25 xl0 -s km/day. 
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100 200 300 1year400days 

Figure 6. Approximate variation of the semi-major axis of Aryabhata 
as a function of time after the launch. 

Hence the above decay in a was expected in the order of magnitude by theoretical 
estimation also. 

(ii) The inclination i has almost remained constant (table 1) throughout the year 
with superimposed periodic variations of low amplitudes. 

(iii) The eccentricity e and the argument of perigee w show large periodic varia¬ 
tions. This may be due to the effect of asphericity of the earth, mainly due to the 
odd zonal harmonics. 

(iv) The nodal regression rate (in right ascension of ascending node Cl) is found 
to be almost constant, having a value of the order of 4-632° per day as expected from 
theoretical calculations. 


Conclusion 

Orbit predictions of Aryabhata were carried out successfully and they met the main 
objective of the mission. Further improvements in the programs developed are 
in progress. When the orbit determination program becomes operational it is 
hoped to carry out a detailed orbit analysis of Aryabhata with the available tracking 
data. However, the orbital information generated was useful in calculating the 
solar-array voltages, temperature conditions of the spacecraft etc. 


Acknowledgements 

We are grateful to Prof. U R Rao and Dr A S Prakasarao for their encouragement 
and fruitful suggestions throughout the period of this work. 




Orbit computations 


301 


List of symbols 


A 

A z 

C D 

CF 

Q,2> $ 2,2 

G 


h 

brain 

Jn 

K v K 2 , K 3 , Ki 


m 

Pn( sin<£), ) 

P 2,2 ( s in 4 1 ) ) 

Re 

r 

r 

s 

U 

U t 

v 

x,y,z 


x,y, z 


x,y, z 


Z 

e 

o a 

A, A, 


F 

p 


Ph 

P T> 

Px j Py> Pz 


Pn> Pei Pz 


<t>> 4>g 


projected area of the satellite 
azimuth 

drag coefficient of the satellite 
control function for visibility 

second order tesseral harmonic coefficients of the earth 

station vector with reference to geocentric inertial co-ordinate 

system 

elevation of the satellite above the ground station 
minimum elevation of the antenna (for Aryabhata h mm = 0) 
zonal harmonic coefficients of the earth of nth order 
factors describing change in atmospheric density due to respectively 
(1) solar radiation intensity F, relative to some average level 
intensity F 0 ; (2) diurnal effect; (3) semi-annual effect and (4) 
geomagnetic perturbation 
mass of the satellite 

Legendre polynomials 

mean equatorial radius of the earth 

radial displacement of the satellite from geocentre 

radius vector of the satellite geocenrtic inertial co-ordinate system 

ballistic parameter of the satellite (s — C D A/2m ) 

geopotential of the earth 

disturbing potential due to asphericity of the earth 
velocity of the satellite 

components of the radius vector with respect to the geocentric co¬ 
ordinate system 

components of the velocity vector in geocentric rotating co-ordi¬ 
nate system 

components of the acceleration vector in the geocentric co-ordinate 
system 

unit vector along the zenith of the station 
local sidereal angle of the station (0 = 6 g + A,) 
sidereal angle of the Greenwich meridian 
geographic longitude of the satellite and station 
gravitational constant of the earth 
atmospheric density 

night time vertical profile of the atmospheric density 
range 

components of the range vector in geocentric inertial co-ordinate 
system 

components of the range vector in topocentric horizontal co¬ 
ordinate system 

geodetic latitude of the satellite and station respectively 
sidereal rate of the earth 
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